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IX.  ABSTRACT  (MAMunum 

A  two-variabla^approach  to  the  model  reduction  problem  with  Hankel-norm  criterion 
is  discussed^  The  problem  is  proved  to  be  reducible  to  obtain  a  two-variable  all¬ 
pass  rational  function,  interpolating  a  set  of  parametric  values  at  specified 
points  inside  the  unit  circle.  A  polynomial  formulation  and  the  properties  of 
the  optimal  Hankel  norm  approximations  are  then  shown  to  result  directly  from  the 
general  form  of  the  solution  of  the  interpolation  problem  considered.  As  a 
consequence,  the  recursive  Nevanlinna  algorithm  can  be  employed  and  the  essential 
stability  properties  of  the  solution  can  be  established  with  the  help  of  the 
Nevanlinna  matrix  (9).  This  short  paper  is  meant  to  briefly  summarize  the  work  in 
the  full  paper  (8),  where  the  reader  is  referred  to  for  mor  details. 
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A  two-variable  approach  to  the  model  reduction 
Procter  with  Hankel  norm  criterion  is  discussed.  The 
prc:!e-  is  proved  to  be  reducible  to  obtain  a  two- 
variable  all-pass  rational  function,  interpolat ing  a  set 
of  para-etric  values  at  specified  points  inside  the  unit 
circle.  A  polynomial  formulation  and  the  properties  of 
the  optimal  Hankel  norm  approximat ions  are  then  shown 
to  result  directly  from  the  general  form  of  the  solution 
of  tre  interpolation  problem  considered.  As  a  conse¬ 
quence,  the  recursive  Nevanlinna  algorithm  can  be  em- 
plo/ed  and  the  essential  stability  properties  of  the 
solution  can  be  established  with  the  help  of  the 
Neva- 1  Inna  matrix  [9].  This  short  paper  is  meant  to 
brief!/  summarize  the  work  In  the  full  paper  [8J ,  where 
the  reader  is  referred  to  for  more  details. 

In ; rcduct ion 

In  1971,  Adamjan,  Arov  and  Krein  ( I ]  have  obtained 
various  powerful  results  related  to  the  theory  of  bounded 
Hanx-I  operators,  which  have  great  significance  in  the 
model  reduction  problem.  Although  their  somewhat  ab¬ 
stract  original  fornulat ion  in  the  framework  of  operator 
trep.-y  has  been  hindering  their  fast  d isseminat ion  into 
t-e  engineering  corn-unity,  the  importance  of  thiseontri- 
buton  ras  begun  to  be  recognized.  The  relevance  of  U] 
tc  :-e  mode!  reduction  problem  was  first  mentioned  by 
R-*5  (7)  in  1572,  while  some  numerical  aspects  of  the 
question  were  published  in  [Jl.  In  (<•],  [5] ,  connect  ions 
be:  -ttr  the  minimal  Mankel-norm  approximat ions  and  the 
o*l*-:ec  state-space  realizations  introduced  by  Hoore 
are  c.t  into  light  and  shown  to  lead  to  an  optimal  ep- 
pre/  -at  ion  algorithm,  requiring  the  solution  of  two 
l-**t-'3/  equations  and  a  singular  value  decomposition. 

f.ung  took  a  pure  one-variable  polynomial  t 
ac  . -cat-  and  proved  that  the  optimal  approximation  prob- 
ic  ti-  so  reduces  to  a  simple  generalized  eigenvalue 
*.:cn  wnen  tne  transfer  function  is  available. 

t  —  r - a’rt-  -o-  ?<t;enal  Aporoxirat len 

'■s  prcslem  ef  rational  approximat ions  In  Hankel- 
-c--  Cw-i  b*  defined  as  follows.  Let  a  t rans far  funet ion 
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h(z)  or  equivalently  Its  Impulse  response  be  given  and 
let  g(z)  bo  some  rational  approximation  of  h(z).  Since 
the  coefficients  of  the  impulse  response  of  a  stable 
linear  system  are  the  Markov  parameters  of  Its  transfer 
function,  comparing  two  Impulse  responses  amounts  to 
comparing  two  Hankel  matrices;  the  Hankel  norm  is  pre¬ 
cisely  the  spectral  norm  of  the  difference  of  these  two 
matrices  and  it  has  a  close  relationship  with  other 
conventional  norms  [6].  The  problem  of  the  optinal 
approximat Ions  can  then  be  stated  In  two  different  forms 
[61. [7], [8]: 

a)  for  a  preassigned  tolerance,  how  to  obtain  a 
minimal  order  approximation  of  h(t) 

b)  for  a  given  order,  how  to  obtain  the  minimal 
norm  approximation  of  h(z) 

Clearly,  the  above  two  problems  are  intimately  related. 

It  turns  out  that  the  problem  of  optimal  approxi¬ 
mations  in  Hankel  norm  can  be  globally  approached  and 
solved  via  a  tw>-variable  polynomial  formalism  [81, 

A.  Two-Variable  Polynomial  Formulation 

let  -  h(z)  -  b(z)/a(z)  be  the  original  transfer 
function  of  order  m. 

-  l(z)  -  zB7(l/z) 

-  r(X,z)t  ir(X,z)  b«  t^>  t*o-v*rIabl«  polyno¬ 
mials  of  degree  n  In  z. 

-  ?(x,z)  »  z'Tu.i/T) 

-  f  can  be  any  constant  of  unit  modulus, 
and  let  us  consider  the  following  equation: 

Xb(z)r(X,z)  -  t  i(z)7(X,z)  •  x(X,z)a(z)  (1) 

The  above  equation  clearly  definas  a  linear  system  of 
equations  in  the  unknown  coefficients  of  the  polyno¬ 
mials  r(\,z)  and  *(X,z).  It  can  be  proved  [8]  that  its 
minimal  degree  solution  r(X,z),  e(X,z)  enjoys  the  fol¬ 
lowing  properties.  For  any  fixed  real  value  of  X  and 
with  (x  (X,z)/Xr(X,z))_  standing  for  tha  stable  projec- 
tion  of  *(X,z)/Xr(X,z)  (l.e.,  after  deleting  its  non- 
strictly  stable  part). 

a)  (e(\ ,z)/Xr (X.r) 1.  is  a  minimal  order  approxi¬ 
mation  o^  hit),  whose  Hankel  norm  does  not 

exceed  IX1*1.' 

b)  If  X*[,  Is  the  (sel)th  singular  value  of  the 

Nankel  matrix  associated  wijh  h(zl  and  If  one 
has  the  strict  ineaual i tv  X  *  X  .  , 
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oisi  ■ t  i«»ao  iffe 
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The  above  results  can  be  easily  proved  with  the 
nelp  of  standard  ratrix  and  polynomial  techniques  by  a 
-ere  reformulation  of  the  problem  [8].  Let  us  indeed 
rewrite  (I)  in  the  following  form: 


i  b(x)  _  r(>  ,z)  r().,i)a (z) 

A  atrr  '  r(>.  ,z)  ”  r(X,z)aiz7 


Let  us  assume  the  zeros  of  a(z)  to  be  distinct  (z-.Zj, 

. ..,z„)  for  the  sake  of  brevity.  With  the  notation 
u.  »  b(z.)/a(z;l,  solving  (I)  is  clearly  equivalent  to 
obtain  tAe  minimal  degree  function  s(X,z)”C?(X,z)/r(\,z) 
interpolat ing  the  parametric  values  Xu.  at  the  points 
a.,  I.e.,  s(X,z.)  -  Xu..  As  a  natter  of  fact,  the 
latter  problem  can  be  solved  [l],(S]  by  the  so-called 
Nevanlinna  algorithm  [s],  which  solves  the  interpolat ion 
problem  In  a  recursive  way.  As  a  result,  the  two-vari¬ 
able  polynomial  r(X,z)  can  be  constructed  recursively 
and  its  stability  properties  investigated  with  the  help 
of  the  Nevanlinna  matrix,  classically  associated  with 
the  Nevanlinna  algorithm  I9j .  Fi»ally,  the  proof  of  the 
minimal  order  and  minimal  Hankel-norm  properties  of 
[-(X,z)/Xr(A.z)]_  is  completed  from  a  direct  link  ex¬ 
isting  between  the  Na/anlinna  ratrix  and  the  singular 
values  of  the  Hankel  -atrix  associated  with  h(z).  [ll, 

[3]. 
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